The binary progenitor of Tycho Brahe's 1572 supernova 
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The brightness of type la supernovae, and their homogeneity as a class, 
makes them powerful tools in cosmology, yet little is known about the 
progenitor systems of these explosions. They are thought to arise when a white 
dwarf accretes matter from a companion star, is compressed and undergoes a 
thermonuclear explosion 1 3 . Unless the companion star is another white dwarf 
(in which case it should be destroyed by the mass-transfer process itself), it 
should survive and show distinguishing properties. Tycho's supernova 4 5 is one 
of the only two type la supernovae observed in our Galaxy, and so provides 
an opportunity to address observationally the identification of the surviving 
companion. Here we report a survey of the central region of its remnant, around 
the position of the explosion, which excludes red giants as the mass donor of 
the exploding white dwarf. We found a type G0-G2 star, similar to our Sun 
in surface temperature and luminosity (but lower surface gravity), moving at 
more than three times the mean velocity of the stars at that distance, which 
appears to be the surviving companion of the supernova. 

Tycho Brahe's supernova (that is, SN 1572) is one of the only two supernovae observed 
in our Galaxy that are thought to have been of type la (the other having been SN 1006) as 
revealed by the light curve, radio emission and X-ray spectra 4-7 . 

The field that contained Tycho's supernova, relatively devoid of background stars, is 
favourable for searching for any surviving companion. With a Galactic latitude b = +1.4°, 
Tycho's supernova lies 59-78 pc above the Galactic plane. The stars in that direction 
show a consistent pattern of radial velocities with a mean value of —30 km s -1 at 3 kpc. 
The predictions of how the companion star would look after the impact, if there is any 
companion, depend on what the star actually is. The star could be in any evolutionary 
stage before the explosion: main sequence, subgiant or red giant 1-3 . The most salient 
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feature of the surviving companion star should be peculiar velocities with respect to the 
average motion of the other stars at the same location in the Galaxy (mainly due to 
disruption of the binary) 8 , detectable through radial- velocity measurements, and perhaps 
also signs of the impact of the supernova ejecta. The latter can be twofold. First, mass 
should have been stripped from the companion and thermal energy injected into it, possibly 
leading to expansion of the stellar envelope that would make the star have a lower surface 
gravity. Second, depending on the interaction with the ejected material, the surface of the 
star could be contaminated by the slowest-moving ejecta (made of Fe and Ni isotopes). 
If the companion's stellar envelope is radiative, such a contamination could be detectable 
through abundance measurements. Therefore, the observations have been designed along 
these lines. The star most likely to have been the mass donor of SN 1572 has to show a 
multiple coincidence: being at the distance of SN 1572, it has to show an unusual radial 
velocity in comparison to the stars at the same location (much above the velocity dispersion 
for its spectral type), and have stellar parameters consistent with being struck by the SN 
explosion. It should also lie near the remnant centre (that is, within our search radius). 

The distance to SN 1572 inferred from the expansion of the radio shell and by other 
methods lies around 3 kpc (2.83 ± 0.79 kpc) 9 . Such a distance, and the light-curve shape 
of SN 1572, are consistent with it being a normal type la supernova in luminosity, as those 
commonly found in cosmological searches 9 . Given the age of the supernova remnant (SNR; 
just 432 yr) and the lower limit to its distance, any possible companion, even if it moved at 
a speed of 300 km s" 1 , could not be farther than 0.15 arcmin (9.1 arcsec) from its position 
at the time of the explosion 8 ' 10 . But the search radius significantly expands owing to the 
uncertainty in the derived centre of the SNR (see Fig. 1). 

We have analysed the stars within a circle of 0.65 arcmin radius, centred on the 
Chandra X-ray Observatory coordinates for the centre of the SNR, up to an apparent visual 
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magnitude V = 22 (Figure 1, Figure 2, Table 1, and Supplementary Tables 1-3). 

All but one of the stars found are either main-sequence stars (luminosity class V) with 
spectral types A4-K3 or giant stars (luminosity class III) with spectral types G0-K3. 

Red-giant stars are possible companions of type la supernovae. Masses in the range 0.9-1.5 
solar masses (0.9-1.5 M ) would be the most favourable cases 11 . Red giants are well 
represented in the sample, but none of them passes the tests for being a viable candidate. 
They are at distances incompatible with that of the supernova. The only giant relatively 
close to the distance of SN 1572 is Tycho A (Fig. 3), but it is closer than SN 1572 and 
shows no peculiarities in velocity spectral type, or metallicity. Main-sequence stars are also 
viable companions of type la supernovae. Close binaries with 2 to 3.5 M Q main-sequence 
or subgiant companions have indeed been suggested as one class of systems able to produce 
type la supernovae 12 . Among systems containing a main-sequence star, recurrent novae 
have been pointed out as possible progenitors 13 . Stripping of mass from the impact of 
the ejecta on this type of companion is also expected 8 ' 14 . Another consequence of the 
impact should be to puff up the star and dramatically increase its luminosity. The size 
and luminosity would later return to their equilibrium values for a star with the new 
decreased mass 8 ' 14 ' 15 . Peculiar velocities should be highest (200-300 km s _1 ) in the case 
of main-sequence companions (orbital separations at the time of explosion are shortest), 
but the measured values of radial velocity (v r ) for the main-sequence stars observed are 
not particularly high. The surface abundances are compatible with solar values. Other 
main-sequence stars (see Fig. 1, Table 1 and Supplementary Table 3) are found at 
wider separations from the geometrical centre, but they have v r values within the range 
corresponding to their respective distances (see Supplementary Discussion). 

We have found a subgiant star ('Tycho G') with lower surface gravity than that of 
main-sequence stars but higher surface gravity than that of red giants, which moves fast in 
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comparison to the mean radial velocities of stars around it, and fits well the expectations 
for distance, reddening and velocity. Comparison of the Tycho G spectrum covering a wide 
wavelength range (3,180-9,400 A) with templates 16 , after dereddening by E(B — V)pz 0.6 
mag, gives a best fit for an effective temperature T e // = 5750 K, a surface gravity log 
g between 4.0 and 3.0, and solar metallicity, which is confirmed by model fitting to 
high-resolution spectra in selected wavelength ranges (see Fig. 3, Supplementary Fig. 1). 
For the spectral type found (G0-G2) and being a slightly evolved star (surface gravity not 
much below the main-sequence value), the mass should be about solar (M 1 M ) and 
thus the radius, for the range of surface gravities above, should be R ps 1-3 R , which 
translates (via our photometric data) into a distance d ~ 2.5-4.0 kpc. This companion could 
have been a main-sequence star or a subgiant before the explosion. While main-sequence 
companions might no longer look like ordinary main-sequence stars after the explosion of 
the type la supernova (and they might resemble subgiants, their envelopes having expanded 
after the supernova impact), subgiants would remain subgiants of lower surface gravity 

9,10,14,15 

Stars at distances d ~ 2-4 kpc, in that direction, are moving at average radial 
velocity 17 v r —20 to —40 km s _1 (in the Local Standard of Rest), with a ~ 20 km s _1 
velocity dispersion 18 ' 19 . Tycho G moves at —108 ± 6 km s -1 (heliocentric) in the radial 
direction. The deviation of Tycho G from the average thus exceeds by a factor of 3 the 
velocity dispersion of its stellar type. It has a 0.3% probability of having that characteristic 
and being unrelated to the explosion (that is, it is a 3a outlier). In contrast, all other 
stars with distances compatible with that of SN 1572 have radial velocities within the 
velocity dispersion as compared with the average of all stars at the same location in the 
Galaxy. We studied through detailed proper motion measurements on the HST WFPC2 
images 20 whether Tycho G has a high tangential velocity as well (see Supplementary Table 
2 and Supplementary Methods). Tycho G has significant proper motion toward lower 
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Galactic latitude: fx b = —6.11 ± 1.34 mas yr _1 (the proper motion along longitude is small, 
Hi = —2.6 ± 1.34 mas yr -1 ). The proper motion in Galactic latitude implies that this star 
is an outlier in proper motion as well, with a derived tangential velocity of 94 ± 27 km s _1 
(a 24 km s" 1 systematic error was added, resulting from a 1.7 mas yr -1 uncertainty in the 
reference frame solution of the images). The other stars do not show such coincidence in 
distance and high tangential velocity. The modulus of the velocity vector has a value of 136 
km s _1 , which is over a factor of 3 larger than the mean velocity value at 3 kpc. 

If Tycho G is the companion star as suggested by its kinematics, the explosion centre 
should have been 2.6 arcsec north of the current location of this star on the basis of its 
velocity. The peculiar velocity would correspond to the peculiar velocities expected from 
the disruption of a white dwarf plus subgiant /main-sequence system 9 ' 10 of roughly a solar 
mass. The system would have resembled the recurrent nova U Scorpii (see Supplementary 
Note 2). The excess velocity corresponds to a period of about 2-7 days, for a system made 
of a white dwarf close to the Chandrasekhar mass plus a companion of roughly a solar mass 
at the moment of the explosion. 

Several paths lead to this star as the likely donor star of SN 1572: its high peculiar 
velocity (both radial and tangential velocities), the distance in the range of SN 1572, and 
its type, which fits the post-explosion profile of a type la supernova companion, as the 
position of this star in the Hertzsprung-Russell diagram is also untypical for a standard 
subgiant. The lower limit to the metallicity obtained from the spectral fits is [M/H] > —0.5 
(see Fig. 3 and Supplementary Fig. 1), which excludes its belonging to the Galactic halo 
population as an alternative explanation of its high velocity. Spectra taken at five different 
epochs also exclude its being a single-lined spectroscopic binary. If our candidate is the 
companion star, its overall characteristics imply that the supernova explosion would affect 
the companion mainly through the kinematics. Our search for the binary companion of 
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Tycho's supernova has excluded giant stars. It has also shown the absence of blue or highly 
luminous objects as post-explosion companion stars. A star very similar to the Sun but 
of a slightly more evolved type is here suggested as the likely mass donor that triggered 
the explosion of SN 1572. That would connect the explosion to the family of cataclysmic 
variables. 
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FIGURE CAPTIONS 

Figure 1 Positions and proper motions of stars. Positions are compared with three centres: 
the Chandra (Ch) and ROSAT (RO) geometrical centres of the X-ray emission, and that 
of the radio emission (Ra). Dashed lines indicate circles of 0.5 arcmin around those centres. 
The supernov position reconstructed from Tycho Brahe's measurements (Ty) is also shown, 
though merely for its historical interest 21 . The radius of the remnant is about 4 arcmin and 
the SNR is quite spherically symmetric, with a fairly good coincidence between radio and 
continuum X-ray emission 8 ' 22 ' 23 . However, there is a 0.56 arcmin displacement along the 
east-west axis between the radio emission and the high-energy continuum in the 4.5-5.8 
keV band observed by XMM-Newton in the position of the western rim 23 (Supplementary 
Note 1). Such asymmetry amounts to a 14% offset along the east-west axis. In SNRs 
from core-collapse supernovae (type II supernovae), up to a 15% discrepancy between the 
location of the compact object and the geometric center is found in the most symmetric 
cases 24 . On the basis of the above considerations, in our search we cover 15% of the 
innermost radius (0.65 arcmin) of the SNR around the Chandra centre of SN 1572. The 
companion star, if there is any, is unlikely to be outside this area (solid line). The proper 
motions of the stars measured from HST WFPC2 images are represented by arrows, their 
lengths indicate the total displacements between AD 1572 and present. Error bars are 
shown by parallel segments. Red circles are the extrapolated positions of the stars back 
to AD 1572. Star Tycho G displays a high proper motion, corresponding to the highest 
tangential velocity in the field, since both stars U and O are at much shorter distances (see 
Supplementary Methods). 

Figure 2 The SN 1572 field and radial velocity of the stars, a, Image from the Auxiliary 
Port at the William Herschel Telescope. It confirms the relative emptiness of the field. 
The search area (see also Fig. 1 bold circle) covers a radius of 0.65 arcmin around RA 
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= 00 h 25 min 19.9 s, dec. = 64° 08' 18.2" (J2000) (the Chandra geometrical centre of 
X-ray emission) with repeated photometric and spectroscopic observations of the included 
stars at various epochs to check for variability and exclude binarity. Additional stars have 
been observed outside of the 0.65 arcmin radius area and are visible in this field (whose 
diameter is 1.8 arcmin). For a remnant distance d = 3.0 kpc and a visual extinction 
Ay = 1.7-2.0 mag toward the candidate stars (see Supplementary Table 3), our search limit 
down to an apparent visual magnitude V = 22 implies that the survey must have detected 
all main-sequence stars of spectral types earlier than K6, plus all subgiant, giant and 
supergiant stars within the corresponding cone. At that distance and with such extinction, 
the Sun would shine as a V = 18.9 mag star, b, Radial velocity (in the Local Standard 
of Rest, LSR) versus distance for the subsample of stars closer than 6.5 kpc (the other 
stars are at a distance well beyond the SNR). We are looking outward along the Galactic 
plane, and the dashed line shows the approximate relationship for the stars in the direction 
of Tycho given by the expression v r = —v Q cos(/ — l & ) + A r sin(2/), where I and /© are 
the respective Galactic longitudes of Tycho and the solar apex, v Q is the Sun's velocity in 
the LSR, and A is Oort's constant 18 . We include two field stars (stars O and U) that are 
slightly away from the search area (at >15% of the radius of the SNR) but at distances in 
the range 2-4 kpc as well. (Star names are labelled lower case in a for clarity.) 

Figure 3 Model fits to observed spectra. Model atmosphere parameters are those listed in 
Table 1, and chemical abundances are solar. They are shown here for our candidate star for 
the companion of SN 1572 (Tycho G) and the red giant (Tycho A) and main-sequence star 
(Tycho B) nearest to the distance of SN 1572 and to the X-ray centre. Identifications of 
the most significant metal lines are given. We have not detected significant spectroscopic 
anomalies, either here or in the whole sample, and most spectra are well reproduced 
assuming solar abundances 25 . Thin lines correspond to the observations and thicker lines 
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to the synthetic spectra. Spectra were obtained at the William Herschel Telescope (WHT) 
with UES and ISIS. Tycho A (bottom panel) is the closest red giant in the sample. It is a 
KO III star, and its mass should be typically M 3 M Q (M Q stands for the mass of the 
Sun). Here, Tycho A is ruled out simply on the basis of having too short a distance. All 
the other red giants are located well beyond Tycho's remnant, and therefore also ruled out 
(see Supplementary Discussion). The A8/A9 star Tycho B (second panel from bottom) 
has M ~ 1.5 M Q , which would fall within the appropriate range for main-sequence type 
la supernova companions, as it would have been massive enough to transfer the required 
amount of mass to the WD. The entirely normal atmospheric parameters, however, strongly 
argue against any such event in the star's recent past. The low radial velocity reinforces 
this conclusion. Tycho G (three upper panels). The second and third spectra from the top 
show computed spectra compared with observed spectra obtained at the WHT with ISIS. 
The upper panel shows the observed spectrum near Ha. This line is blueshifted, implying a 
peculiar radial velocity exceeding about 3 times the velocity dispersion for its stellar type. 
This star does not belong to the halo population (Supplementary Fig. 1). 
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TABLE CAPTION 



Supernova companion candidates within the search radius and limiting magnitude. Angular 
distances 6 are from the Chandra X-ray geometrical centre, located at RA = 00 h 25 
min 19.9 s, dec. = 64° 08' 18.2" (J2000). Synthetic spectra, under the assumption of 
local thermodynamic equilibrium (LTE), are fitted to the observed ones using the grids of 
model atmospheres and the atomic data of Kurucz 25 , with the Uppsala Synthetic Spectrum 
Package 27 . This determines the atmospheric parameters effective temperature T e ff and 
surface gravity g. Intrinsic colours and absolute visual magnitudes are deduced from the 
relationships between spectral type and colour and between spectral type and absolute 
magnitude for the different luminosity classes 28 . Comparison with our photometric BVR 
measurements (see Supplementary Table 3) yields the reddening E(B — V), from which 
the visual extinction Ay and the corrected apparent visual magnitude V are calculated. 
Comparison with the absolute visual magnitude then gives the distance d. Uncertainties 
in T e ff are 250 K. Tycho J is a binary of main-sequence stars with masses in the range 
0.80-0.85 M and quite similar atmospheric parameters. Tycho C is found in HST images 
as being two stars (CI and C2) 0.25 arcsec apart. Modelling of the composite spectrum 
and the HST magnitudes of the stars show that they do not constitute a physical binary; 
the hot fainter component C2 is at larger distance than CI. Within this list of stars, D, G, 
N and V have proper motions along Galactic longitude and latitude of Hi = -3.23, ^ = 
-0.58 ± 0.66 (same error in both coordinates, units in mas yr _1 ) for star D, /// = -2.60, fib 
= -6.11 ± 1.34 for star G, m = 3.23, fi b = 1.45 ± 1.15 for star N, and ^ = 1.61, jj b = -2.85 
±0.78 for star V. 
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Table 1: Table 1. Characteristics of the supernova companion candidates 



Star 


e 

(arcsec) 


Spec, type 
& lum. class 


T e// 

(K) 


log g 
(c.g.s.) 


E(B-V) 
(mag.) 


d 

(kpc) 


Tycho A 


1.6 


K0-K1 III 


4750 


o c+0.5 
z -°-0.5 


U.OO_o.05 


1 1+0.3 
i - i -0.3 


Tycho B 


1.5 


A8-A9 V 


7500 


4 -°-0.5 


o.6ol8:8i 


2.618:1 


Tycho CI 
Tycho C2 


6.5 
6.5 


K7 V 
F9 III 


4000 
6000 


4.518:1 

2.018:1 


o.5l8:l 
o.6l8:l 


o.75l8:l 

> 20 


Tycho D 


8.4 


Ml V 
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4 c;+o.5 

4 -°-0.5 


0.618:1 


o.8l8:l 


Tycho E 
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2.ol8:l 
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Tycho F 


22.2 


F9 III 
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o.54i8:l 
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Tycho G 


29.7 
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q pr+0.5 
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o n+1.0 
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Tycho H 


30.0 
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5000 


3.ol8:l 
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> 13 


Tycho J 


33.9 


Kl V 


5000 


4-518:1 


0.5818:}? 


-x . _i_n s 
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Tycho K 


35.0 
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6000 


2.0181 
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> 10 


Tycho N 


35.4 


GO V 


6000 


4.5181 


o.62i8:8? 


o 1+0.7 
z - i - 0.7 


Tycho V 


29.2 


K3 V 
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4.5181 


o.6ol8:18 


q O+0.6 

<J-O-0.6 
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Fig. 1.- 
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Supplementary Methods 

Specification of the WFPC2 measurements 

The proper motions were calculated by comparing through a maximum likelihood 
analysis the displacement in image centroids of the target stars to a reference frame 
determined from all the point sources in common between the WFPC2 images of the two 
epochs 1 . 

For the majority of our targets, the stellar images fell on the WF3 CCD of the WFPC2 
mosaic camera on both epochs. However, due to differences in the camera pointing and 
roll-angle between the exposures at the two epochs, the target star K fell on different CCDs 
of the WFPC2 mosaic at the two epochs, as did the stars L, M, and W (these three not 
being targets). The overlap region proved to be too small to determine a reference frame, 
and so the proper motions of these stars could not be measured. A similar effect entered 
the measurement of target H, which had to be done using images spanning a shorter 
time baseline, and making the uncertainty too large. No accurate image centroids were 
derived for the stars CI and C2 due to their small separation, neither for Target A which is 
saturated in the reference images and has a background star at half a second of its centroid. 
The other stars do not have significant proper motions (at > 2a). This proper motion 
programme continues in HST Cycle 13 where measurements with smaller error bars will be 
obtained using both WFPC2 and ACS. 

1 Ibata, R. A. & Lewis, G. F., Proper motion measurements with WFPC. Astron. J, 116, 
2569 (199 
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Supplementary Table 1. Spectroscopic runs in this work 



Date 




Telescope 


Instrument / Grating 


Resolution 


Coverage (A) 


1997- 


■12-04 


WHT 


ISIS+R158B/R158R 


600 


Qonn 
ozUU- 


nnnn 

-yuuu 


1998- 


■01-03 


WHT 


ISIS+R300B/R158R 


1700-600 

X 1 WW www 


oonn 


nnnn 

yuuu 


1999- 


•02-24 


WHT 

V V -L J. -L 


ISIS+R300B/R158R 


1 700-600 

-L 1 WW www 


oZUU- 


nnnn 
-yuuu 


2000- 


•11-21 


Keck I 


ESI 


10000-7000 


/innn 

4UUU~ 


i nnnn 

1UUUU 


2001- 


•06-16 


WHT 


UES 


50000 


a nnn 

4UUU- 


71 nn 

- ( 1UU 


2001- 


•06-17 


WHT 


ISIS+H2400B /R1200R 


10000-7000 


/innn 
4uuu- 


7nnn 
- ( uuu 


2002- 


•06-12 


WHT 


ISIS+H2400B /R1200R 


10000-7000 


44nn 


7nnn 
- ( uuu 


2002- 


•06-13 


WHT 


ISIS+ H2400B /R1200R 


10000-7000 


a /inn 


7nnn 
- ( uuu 


2002- 


•06-14 


WHT 


ISIS+R300B/R158R 


1700-600 


ozuu- 


nnnn 
-yuuu 


2002- 


•10-09 


NOT 


ALFO S C + Grism9 


4500-3900 


/i/inn 


annn 
-yuuu 


2002- 


•12-05 


WHT 


ISIS+H2400B /R1200R 


10000-7000 


a cnn 


vnnn 
/UUU 


2002- 


•12-06 


WHT 


ISIS+H2400B/R1200R 


10000-7000 

-LWWWW 1 WWW 


4400- 


-7000 


2003- 


•07-27 


Keck II 


LRIS+300/5000 


850 


3200- 


-7900 


2003- 


•07-30 


WHT 


ISIS+R1200B/R1200R 


7000 


4000- 


-7000 


2003- 


•07-31 


WHT 


ISIS+R1200B/R1200R 


7000 


4000- 


-7000 


2003- 


-11-16 


WHT 


ISIS+R1200B/R1200R 


7000 


4000- 


-7000 


2003- 


■11-17 


WHT 


ISIS+R1200B/R1200R 


7000 


4000- 


-7000 


2003- 


■11-18 


WHT 


ISIS+R1200B/R1200R 


7000 


4000- 


-7000 


2003- 


•12-20 


Keck I 


LRIS+150/7500 


500 


3156- 


-9400 



Spectroscopic runs at the 4.2 m William Herschel Telescope (WHT), the 2.5 m Nordic Optical 
Telescope (NOT) and the 10 m Keck I and Keck II telescopes used in this work 
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Supplementary Table 2. Hubble Space Telescope data sets used for proper motion 

measurements 



Dataset 


Date 






Filter 


Exposure Time 


Instrument 


U8S9010HM 


2003- 


11- 


•08 


F555W 


100.0 s 


WFPC2 


U8S9010IM 


2003- 


11- 


•08 


F555W 


100.0 s 


WFPC2 


U8S9010JM 


2003- 


11- 


•08 


F555W 


100.0 s 


WFPC2 


U8S9010KM 


2003- 


11- 


•08 


F555W 


100.0 s 


WFPC2 


U8S9010LM 


2003- 


11- 


-08 


F555W 


100.0 s 


WFPC2 


U8S90201M 


2003- 


11- 


-06 


F555W 


400.0 s 


WFPC2 


U8S90202M 


2003- 


11- 


-06 


F555W 


400.0 s 


WFPC2 


U8S90203M 


2003- 


11- 


-06 


F555W 


400.0 s 


WFPC2 


U8S90204M 


2003- 


11- 


-06 


F555W 


400.0 s 


WFPC2 


U8S90205M 


2003- 


11- 


-06 


F555W 


1.0 s 


WFPC2 


U8S90206M 


2003- 


11- 


-06 


F555W 


1.0 s 


WFPC2 


U8S90207M 


2003- 


11- 


-06 


F555W 


1.0 s 


WFPC2 


U8S90208M 


2003- 


11- 


-06 


F555W 


1.0 s 


WFPC2 


U8S90209M 


2003- 


11- 


-06 


F555W 


1.0 s 


WFPC2 


U8S9020AM 


2003- 


11- 


-06 


F555W 


1.0 s 


WFPC2 


U8S9020BM 


2003- 


11- 


-06 


F555W 


0.2 s 


WFPC2 


U8S9020CM 


2003- 


11- 


-06 


F555W 


0.2 s 


WFPC2 


U8S9020DM 


2003- 


11- 


-06 


F555W 


0.2 s 


WFPC2 


U8S9020EM 


2003- 


11- 


-06 


F555W 


0.2 s 


WFPC2 


U8S9020FM 


2003- 


11- 


-06 


F555W 


0.2 s 


WFPC2 
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Dataset 


JJatc 






r liter 


Exposure Time 


Instrument 


TTocnm nui\ ft 


onno 


11- 


no 
-Uo 


TTir rr rlTT 
r 000 W 


1 nn n ^. 
1UU.U S 


Wrruz 


TTocnm nTA ft 

U8d9U1U1M 


z(J(Jo- 


i i 
11- 


no 

-08 


TT 1 rr rr rr T * 7 

b 555 W 


1 nn n ^ 
1UU.U S 


Wr rU2 


TTocnm n tat 

U8b9UlUJM 


anno 

z(J(Jo- 


11- 


no 

-08 


TTi rr rr rr "\ s 7" 

1* 555 W 


1 nn n ^ 
1UU.U S 


Tl 7T7 1 T-)00 

Wr rU2 


TTocnm riT/A/T 
UooyUlUKM 


onno 
ZllUo- 


11- 


no 
-Uo 


tpf: rr cats/" 
r 555 W 


1 nn n ^. 
1UU.U S 


AS7TT1T300 


TTOCnm TTT AT 

U8b9UlULM 


onno 

z(J(Jo- 


11- 


no 

-Uo 


r 555 W 


1 nn n ^. 

lUU.U s 


Tl 7T7 1 T-)00 


TTocnnom at 


zUUo- 


11- 


-Uo 


T?rr. rr z.~\xt 
r 555 W 


4UU.U s 


AS7TT1T300 


TTocnnonoi\ ft 

U8b9U2U2M 


anno 

z(J(Jo- 


11- 


-Uo 


T7 1 rr rr rr T * 7 

r 555 W 


4UU.U s 


Tl 7T7 1 T-)00 


TTocnnonoi\ ft 

U8b9U2U3M 


onno 

z(J(Jo- 


11- 


-Uo 


T7 1 CCCHT 

r 555 W 


4UU.U s 


Tl 7T7 1 T-)00 


TTOOnnon/IAT 

U8b9U2U4M 


z(J(Jo- 


11- 


-Uo 


TTirr rr rrT*7 

l 1 555 W 


/|AA A _ 

4UU.U s 


Tl 7T7 1 T-)00 


TTocnnonci\ ft 

U8b9U2U5M 


z(J(Jo- 


■11- 


-Uo 


T7irr rr rrT* 7 

.b 555 W 


1 n ^ 

l.U s 


WrrU 


TTocnnonri\ ft 

U8b9U2U6M 


z(JUo- 


■11- 


np 

-Uo 


TTirr rr rrT*7 

i 1 555 W 


1 n ^ 

l.U s 


WrrU 


t Toonnon7i\ ft 


zUUo- 


11- 


n^ 
-Uo 


TJ'rr rr rr"** 7 

r 555 W 


1 n 

l.U S 


AS7TT1T300 


TTocnnonoi\ ft 

U8b9U2U8M 


2UU3- 


11- 


-Uo 


T7 1 tr rr rr \it r 

1 555 W 


1 n ^ 
l.U S 


WrrU 


TTocnnonni\ ft 


onno 

2UU3- 


11- 


n/^ 

-Uo 


T7 1 rr rr rr 7 

b 555 W 


1 n ^ 

l.U s 


WrrU 


U8b9U2UAM 


2UU3- 


11- 


n£? 

-Uo 


T7 1 tr rr rr \it 7 

r 555 W 


1 n ^ 

l.U s 


AS7T7 1 T)00 

Wr rbz 


TToonnonTDAT 
UobyUzUBM 


onno 
zUUo- 


11- 


n^ 
-Uo 


TPrr rr rrTi 7 

r 555 W 


n ^ 
0.1 S 


AS7TT1TDOO 


TT8C!Qn9nPl\/T 

uooyuzuoivi 


zuuo- 


1 1 
ii- 


-uo 


r 000 vv 


n 9 c 


WPPP9 
VV r Jr 


U8S9020DM 


2003- 


11- 


■06 


F555W 


U.2 s 


WFPC2 


U8S9020EM 


2003- 


11- 


■06 


F555W 


U.2 s 


WFPC2 


U8S9020FM 


2003- 


11- 


-06 


F555W 


U.2 s 


WFPC2 


U8S9020GM 


2003- 


11- 


-06 


F555W 


U.2 s 


WFPC2 


U8S9020HM 


2003- 


11- 


-06 


F555W 


1UU.U s 


WFPC2 


U8S9020IM 


2003- 


11- 


-06 


F555W 


1UU.U s 


WFPC2 


U8S9020JM 


2003- 


11- 


-06 


F555W 


1UU.U s 


WFPC2 


U8S9020KM 


2003- 


11- 


-06 


F555W 


1UU.U s 


WFPC2 


U8S9020LM 


2003- 


11- 


-06 


F555W 


1UU.U s 


WFPC2 



Supplementary Table 2 (cont.) The data from 1999 correspond to HST programmes 
G06435 and G074U5. The data from 2UU3 correspond to G09729. 
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Supplementary Table 3. Magnitudes and radial velocities of the stars 



Star 


9 


B 


V 


R 


v r 




(arcsec) 


(mag.) 


(mag.) 


(mag.) 


(km/s) 


A 


1.6 


14.8218;" 


13.2918:81 


1 9 24 +0 - 03 

lZ.Z+_ 03 


-2311 


B 


1.5 


16.351831 


i5.4H8:81 


15.1118:18 


-381^ 


C* 


6.5 


2i.06i8il 


19.0618:8^ 


17.7718:81 


-3318 


D 


8.4 


22.9718:1 


20.7018:18 


19.3818:88 


— 


E 


10.6 


21.2418:11 


19.7918:8? 


is.84±8:8l 


-2611! 


F 


22.2 


19.0218:^ 


i 7 7 q+0.03 
L ' ■ ' °-0.03 


16.9418:81 


-341H 


G 


29.7 


20.0918:8^ 


18.7118:81 


17.8318:81 


-9918 


H 


30.0 


21.39l8il 


19.8018:8? 


is.78l8:8i 


-7H18 


J 


33.9 


2i.i5l8i| 


19.7418:8? 


18.8418:8^ 


-4518 


K 


35.0 


21-6418:11 


on 1 1 +0.08 
ZU.±±_q og 


i9.i5l8:8i 


-33118 


N 


35.4 


19.5918:88 


18.2918:81 


1 7 47+0.03 

1 ' - 4 ' -0.03 


-30I8 


V 


29.2 


oo oo4-0.33 
ZO.OZ_ .33 


21 41 + - 13 
zi.+i_ 13 


2o.2o+g : gl 


-47118 





41.5 


18.6218:81 


1 7 9 q+0.03 
1 1 .ZO_ 03 


16.3718:81 


-151? 


P 


40.4 


18.8418:18 


17 fi i+0.03 
1 ' - Oi -0.03 


16.7818:81 


-36118 


u 


39.5 


i9.03l8:8t 


1 7 7Q+0-03 
1 ' • ' °-0.03 


16.9518:81 


-3811 



Angular distances from Chandra's geometrical X-ray centre, BVR apparent magnitudes 
and radial velocities (LSR) for the sample of SN companion candidates. Radial velocities 
have been measured from the wavelength shifts of several absorption lines in each observed 
spectrum (see Supplementary Table 1). *Data are for the unresolved pair. From the HST 
data, the brighter, bluer component has magnitudes B = 21.28, V = 19.38, R = 18.10 while 
the fainter, redder component has B = 22.91, V = 20.53, R = 19.23. 
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Supplementary Note 1 

Note on the origin of the asymmetry in Tycho SNR 

Evidence that the eject a encountered a dense H cloud at the eastern edge giving rise 
to brighter emission and lower ejecta velocity there, while finding a lower-density medium 
in the western rim, might account for the asymmetry 1 . 

To obtain an approximate position of the dynamical center of the explosion one 
would need to trace back the expansion of the SNR to the time of the explosion from the 
SNR expanding filaments. Unlike in SNe II, neither in SN 1572 nor in any other SN la 
has the derivation been possible thus far, due to the faintness of the filaments. This 
measurement remains an interesting long-term project to be undertaken with the use of 
large telescopes. In this work, arguments based on what is found in SNe II together with 
the E-W geometrical asymmetry of this SNIa are used to establish the area of the search. 

1 Decourchelle, A. et al. XMM-Newton observation of the Tycho supernova remnant Astron 
& Astrophys. 365, L218-L224 (2001). 
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Supplementary Discussion 

Further discussion on red giants and main sequence stars 

Red giants In the red-giant companion case, the envelope is loosely bound gravitationally, 
and upon collision with the SN la ejecta it should be either completely stripped or just a 
small fraction of it remain bound to the core 1,2 . If the full envelope were lost, the remaining 
He core would appear as a hot He pre- white dwarf, not as a red giant. Our survey would 
have revealed the presence of hot helium stars resulting from complete stripping of the 
envelope of a red giant because they are brighter than the other stars that we do observe. If 
a minor fraction of the envelope of the red giant were retained instead, the H-burning shell 
would remain active and the residual envelope would expand to red-giant size. In the case 
of Tycho A, whose envelope is ~ 2.6 M Q , nothing similar to the star that we are now seeing 
could ever be what remains after the supernova impact. Tycho A and the other red giants 
are in any case ruled out on the basis of their incompatible distance with SN 1572. 

Main sequence stars We should emphasize that there exists, beyond differences in the 
post-explosion evolution calculations, unanimous agreement that the companion object, 
if it were a main-sequence star, should exhibit an odd combination of stellar parameters 
(log g and T e ff) for its mass and that it should move fast 1,3 ' 4,5 . That is not found in the 
observed sample of main-sequence stars. Apart from Tycho B (see Figure 3), there are six 
main-sequence stars in Table 1. One of them (Tycho J) is a spectroscopic binary made of 
two main-sequence stars with similar spectral types (G8/K0-K3). Its radial velocity fits the 
derived distance, and no such system is likely to be the surviving companion of any SN la. 
As stated above, all main-sequence stars of spectral type earlier than K6 were detected to 
the distance of the remnant. The total mass available for transfer in main-sequence stars 
of type later than K6 (M <,0.6 M ) excludes these stars as viable candidates to SN la 
companions. 
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Supplementary Figure 1. Spectra of Tycho G showing it to be a subgiant not belonging 
to the halo population. (Right panel) Several fits to Fe and Ni lines in Tycho G for solar 
abundances (bold) and abundances [Fe/H] = —0.5 (dashed line) and [Fe/H] = —1 (dotted 
line). This star does not belong to the halo population: it shows solar metallicities in 
Fe and Ni. The observed spectra were obtained with ISIS at the WHT. (Left panel) A 
low-resolution spectrum over a wide wavelength range was obtained with LRIS at the Keck 
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Observatory (second from top) and it is compared with template model spectra of the same 
spectral class and various metallicities. It has also been used to determine the atmospheric 
parameters (see main text). Very high and low metallicities are excluded. The spectrum of 
Tycho G confirms that it is a subgiant star. 
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Supplementary Note 2 

Tycho Brahe SN as a U Sco system 

The evolutionary path that gave rise to SN 1572 could be similar to that leading to the 
recurrent nova U Scorpii pointed out as a candidate to SN la progenitor. The companion 
candidate is now a mildy evolved star within the solar mass range. The excess velocity 
implies an orbit of 6 days before the explosion. Starting from a white dwarf with a mass 
~ 0.8 Mq plus a somewhat evolved companion of ~ 2.0-2.5 Mq filling its Roche lobe, and 
with a period of 12 days, it would have ended up as a white dwarf at the Chandrasekhar 
mass (~ 1.4 Mq) plus a companion of roughly 1 Mq, the period then being ~ 6 days 
(orbital velocity ~ 90 km/s). The effective radius of the Roche lobe of the companion just 
before explosion would have been ~ 7 Rq. Now the radius of the companion has to be less 
than three times the solar radius given the effective temperature and luminosity of the star. 
The smaller current radius would result from the opposite effects of mass stripping and 
shock heating by the supernova impact, plus subsequent fast cooling of the outer layers up 
to the present time. 



